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Abstract 

In this article we construct a tiny graviton matrix model for type 
IIB string theory on a plane-wave background with null dilaton. For 
the linear null dilaton case, we analyze its vacuum and the excitation 
spectrum around the vacuum, and discuss the time-dependent fuzzy 
three-sphere solutions and their evolution. It turns out that at very 
late time the non-Abelian fuzzy degrees of freedom disappear, which 
indicates the appearance of perturbative strings. 
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1 Introduction 



One of central problems in the study of string theory is to look for its 
non-perturbative definition. One decade ago, Banks, Fischler, Shenker and 
Susskind (BFSS) suggested that a matrix model could be a candidate for 
non-perturbative definition of string theory (M theory) pQ. They constructed 
so-called BFSS matrix model on the flat background, and conjectured that 
large N limit of their model is equivalent to uncompactified M theory in 
the infinite momentum frame. In (2J [3], the authors showed how to obtain 
BFSS model from discrete light-cone quantization (DLCQ) of M theory. The 
BFSS model is a (0 + l)-dimensional super- Yang-Mills gauge theory, or a 
matrix model, and its Hamiltonian can be obtained from the nonrelativistic 
Hamiltonian of D-particles. Remarkably, there is another way to obtain this 
model by doing matrix regularization of super-membrane in flat space-time. 

Since the construction of BFSS model, people has tried hard to construct 
the matrix models on general backgrounds. For the weakly curved spacetime, 
the authors in [I] studied the non-Abelian D-particle action, which could be 
essential to get the matrix model action. Even for this case, the complete 
D-particle action is still not clear partially due to notorious difficulty in 
constructing non-Abelian D-brane action. However, for one class of curved 
spacetime, so-called plane-wave background, the corresponding BMN matrix 
model was well-established and much better understood [5]. Similar to the 
matrix model in flat spacetime, the BMN model could be obtained from the 
matrix regularization of supermembrane on the background [6]. This suggest 
that the matrix regularization of super-membrane could be an important 
way to construct the matrix model in curved spacetime. For the general 
curved spacetime, the action of the supermembrane could not be determined 
completely [TJ. However, it is remarkable, for a general class of plane- wave-like 
background, the matrix model could be constructed to all order [8]. 

There are also other nonperturbative definition of various kinds of string 
theory. For IIA theory, one non-perturbative definition of string theory is 
matrix string theory [9]. For IIB theory, one interesting construction is the 
IKKT matrix model [TU], which is a (0 + 0)-dimensional theory. Quite re- 
cently, a (0 + l)-dimensional matrix model associated to IIB theory on plane- 
wave background was proposed by Sheikh- JabbaripTjS. His observation is 
that the matrix models, including BFSS and BMN, can be understood as 
the theory of gravitons. In BFSS case, the D-particles are gravity waves in 
11- dimensional flat spacetime, while in BMN case, the fundamental degrees 
of freedom are tiny graviton, which are membrane gravitons with very small 

■'■For an alternative proposal to the Tiny graviton matrix theory, see |12j . 
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size. In other words, the DLCQ of M-theory in 11-dimensional background 
could be described by a theory of N membrane graviton. This led Sheikh- 
Jabbari to conjecture that DLCQ of IIB superstring theory on plane-wave 
background is a + 1-dimensional matrix model, whose Hamiltonian can be 
written down through the action of tiny gravitons. In this case the tiny gravi- 
tons are very small spherical D3-branes, so the Hamiltonian can be obtained 
from light-cone matrix regularization and quantization of D3-brane Hamilto- 
nian on the background. More study about tiny graviton matrix model can 

be found in PI EH H5] . 

The search for a non-perturbative definition of string theory is not only 
just to complete the formulation of the theory, but may also has profound 
physical applications. In particular, at the very early universe, the appro- 
priate description of the physics there could be in terms of nonperturbative 
string degrees of freedom since the string coupling could be strong. One 
concrete realization of this idea is so called matrix big bang model and its 
cousins. The original matrix big bang is starting from the IIA string theory 
in null dilaton background |16j. Formally, this IIA string theory is exactly 
solvable perturbatively but since the string coupling becomes strong at null 
infinity, the perturbative string theory breaks down. It turned out that the 
matrix degrees of freedom, rather than the point particle or the perturbative 
string, could describe the physics near the null singularity. Following the 
prescription in [21 E] , a dual matrix string theory was proposed. The study 
of the plane-wave matrix big bang model via DLCQ was carried out very 
recently [17]. For a nice review on this model, its generalization and relevant 
references, please see [18] . 

In [19] , one class of solvable IIB string theory in linear null dilaton back- 
ground has been investigated. The background is 

ds 2 = —2dx + dx~ — \(x + )x 2 dx + dx + + dx 1 dx 1 , 

(j) = 

(^) + 1234 = (i ? 5 ) + 5678 = 4/(2 + ), (1.1) 

in which one has to require 

A(x + ) = -\<P"{x + ) + f(x+)e 2 ^ + \ (1.2) 

to satisfy the equations of motion. This background keeps sixteen non- 
linearly realized supersymmetries. The background has continuous symmetry 
algebra so(4) © so(4) (B s h(8), in which two so(4) are the residual rotational 
symmetry algebra and h(8) is the Heisenb erg-type algebra. 
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Generically the background (11.11) has curvature singularity. To show this 
clear, it is better to work in the Einstein frame. For the linear null dilaton 
case we are interested in this paper, 

(f)=-cx + } with c > (1.3) 

the issue has been discussed carefully in section two of [19]. It has been 
shown that no matter what the background flux is the background ( 11. ip is 
geodesically incomplete and x + — * — oo is a curvature singularity where an 
inertial observer will feel a divergent tidal force. The existence of singularity 
and the fact that the string coupling is also divergent at x + — > — oo both 
suggest that the perturbative string degrees of freedom to which the graviton 
belong can not give correct description of the physics there, and instead 
nonperturbative string will take over and play a fundamental role. In other 
words, to describe the physics near x + — > — oo calls for a nonperturbative 
definition of string theory, which might be a matrix model. 

Similar to the usual plane-wave background, one can do light-cone quan- 
tization of perturbative string in the background ( 11.11) . However, due to the 
absence of translational Killing vector along x + , the lightcone Hamiltonian is 
not conserved. Actually both the Hamiltonian and perturbative string states 
are time-dependent. Nevertheless, because this background is invariant un- 
der translation in the x~ direction, the compactification of x~ on a circle 
gives a conserved lightcone momentum. This fact also suggests that DLCQ 
of string and D-brane in this background is feasible. In this paper, we will 
follow the proposal in [11] to construct a time-dependent matrix model on 
the background (11.11) . 

This article is organized as follows. In section 2, after constructing the 
matrix model, we analyze the spectrum around the vacuum, and compare 
its spectrum with the spectrum of perturbative string on this background. 
In section 3, we look for and discuss the fuzzy three sphere solutions. In our 
cases, the fuzzy three-sphere solutions are time-dependent, with their radius 
varying with time. At last, we end with conclusions and some discussions. 
In the Appendix, we discuss the giant graviton in the background and show 
that its dynamics is the same as that of fuzzy three sphere. 

2 The matrix model 

To construct the matrix model, we follow the construction proposed in 
|llj . Let us start from the light-cone Hamiltonian of D3-brane. The Born- 
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Infeld and Chern-Simons action for D3-brane ig^| 

S = J dr&oe-* yj '-det(G^d^X»d u ,X») + J C 4 . (2.1) 

Here G^, C4 are the spacetime metric and Remond-Remond four-form po- 
tential of the background (II. ip . The X^'s are the embedding coordinates. 
The indices n',v' = 0,1,2,3 label the D-brane's world- volume coordinates, 
and i',f, ■ ■ ■ denote the spatial ones. 
In the light-cone gauge, 

X + = r, G lu ,d Q X"d v X v = 0. (2.2) 

We use P_ to denote the conjugate momentum of X~, P + = — P_, and Pi 
to denote the conjugate momenta of X 1 , I — 1, 2, . . . , 8. Using the light-cone 
Hamiltonian density 7i = P-d T X~ + PjdrX 1 — C, we have the bosonic part 
of light-cone Hamiltonian: 

H bos = J d 3 a[^(P? + Pi) + ^~(Xf + Xl) 

p -2*(t) I 

+ ^P^ ( 3! {X *> X *> Xk]{X% Xj > Xk} 
+ ^{X a ,X b ,X c }{X a ,X b ,X c } 

+ ^{X\ X'\ X a }{X\ X\ X a } + ^{X\ X\ X b }{X\ X\ X b }) 
) ^Hhix^X 3 , X k , X 1 } + e abcd X a {X b , X c , X d })], (2.3) 
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where = 1, 2, 3, 4, a, b, ■ ■ ■ = 5, 6, 7, 8. In the above relation the Nambu 

three brackets are used, 

{X 1 , X J , X K } = j' j ' k 'a i ,X I d j ,X J dvX K . (2.4) 

The fermionic part of the Hamiltonian can be obtained through super- 
space techniques [2Ul |2"T| |2"2"]. On the background (11.11) . it is given by 

H fer = J d 3 a[fe^ af3 ^ 

+^^ aP (° ij )i{x i , x\ M + ^{o ab t{x\ x b , 



2e' 

+^(a ab )i{X^X b ,^})]. (2.5) 



2 We will always set l s = 1 in this paper. 
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Here the ^'s carry two spinor indices, corresponding to two SO(4) (SU (2) p x 
S77(2)_r)'s Weyl spinor representations respectively. 

To do discrete light-cone quantization, the light-like direction should be 
compactified. 

X~ = X~ + R^. (2.6) 
With this compactification, P + is quantized, 

P + = , (2.7) 

where J G Z. 

Next one has to quantize Nambu three brackets in order to have a ma- 
trix regularization of the D3-brane Hamiltonian. It is well known that the 
quantization of Nambu odd brackets is very difficult, because some impor- 
tant properties of Nambu odd brackets will get lost if we apply the usual 
way of quantization of Nambu even brackets. In [TT], Sheikh- Jabbari pro- 
posed a way to quantize the Nambu odd bracket by introducing a special 
constant matrix. For Nambu three brackets, the constant matrix is £5, and 
the quantization is to replace {F, G, K} with J[Fj x j, Gj x j, Kj x j, £ 5 ], where 

[Fi, F 2 , F 3 , F A ] = ±e iM F h F i2 F i3 F i4 . (2.8) 

The £5 is made out of the direct product of unit matrix J 4X 4 and SO (4) 
chirality operator T 5 . In this paper, we only need one property of it. Namely 
there exist matrices X l JxJ , i = 1,2,3,4, which satisfy 



p2 

r y i yi v k r 1 — Ji kl y 1 



E^xi) 2 = R 2 h,j- (2.9) 

i=l 

Recalling that a three sphere of radius R can be embedded in four dimensional 
space as 

{X\X 3 ,X k } = R 2 e ijkl X l , 
4 

J^pT) 2 = R 2 , (2.10) 

i=i 

we know that the matrices satisfying f!2.9p is actually a realization of fuzzy 
three sphere Sp of radius R. For the explicit construction, please see [TT] 
and [231. 
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H = R-Tt{-U 2 j + --^X] + ^[X 1 , X J , X K , £ 5 ] [X 1 , X J , X K , C 5 



Now for the sector with J units of light-cone momentum, we replace 
X 1 , P 7 , ^ ay 3, ijj^h with J x J matrices. 

ij) v/J^Jxj. (2.11) 

In the following part of this paper, we will denote Xj x j, Hj x j, ipj x j with 
X 1 , II 7 , ^. We wish this will not bring any confusion. 

With all these preparation, we can finish constructing the matrix model 
for DLCQ of IIB string theory on background (11.11) by doing the replace- 
ments: 

J d 3 a -> P„Tr (2.12) 

{F, G, K} —> J[F, G, K, £ 5 ]. (2.13) 
Then the matrix model Hamiltonian is 

1 TT2 1 A V2 e ~ 2 ^ 

-^-(e llkl X l [X\ X\X\ £ 5 ] + e abcc X a [X 6 , X c , X d , £ 5 ]) 

+ Z^W^ _ ^ tA4 ) + 2e^(^ ta/3 (^)i[X\ X', ^, £ 5 ] 

+^(^)^[x a ,X b ,^,£ 5 ]) - 2 e ^(^(^)|[X i ,X^^,£ 5 ] 
+^ ab t\X a , X b C,])} (2.14) 
In this Hamiltonian, A, <p, f depend on r, and satisfy 

A(r) = -^"(r) + / 2 (r)e 2 ^. (2.15) 

Now we have constructed a (0 + 1) -dimensional matrix model of type IIB 
string theory on time-dependent background (11.11) . The symmetry algebra 
of background (11.11) is [so(4) © so(4)] © s ^(8), where h(S) is Heisenberg- 
type algebra. So DLCQ of string theory on this background should respect 
50(4) x 50(4) symmetry. From (12.14)) . it is easy to see that our matrix 
model do have this symmetry. 

Generically the matrix model has time-dependent mass-squared A and the 
coupling constant is also time-dependent. In this paper, we are interested in 
the null linear dilaton background 

(j) = -ct, c> 0. (2.16) 
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so that 

A = /V^. (2.17) 
As the usual PP-wave case, the string theory limit could be 

J, -> cx), with P+ being fixed. (2.18) 

In this limit, our matrix model should describe light-cone type IIB string 
theory on the background (11. ip . 

Because our matrix model is time- dependent, it is difficult to find a static 
solution. The unique static solution is the vacuum solution: 

X 1 = X a =ip a/3 = = 0. (2.19) 

Let us consider the fluctuations around the vacuum. Expanding the 
Hamiltonian to the second order of the fluctuations around the vacuum, 
we obtain 

H% = Tr[(^n t 2 + ^-Xl) + (^Ul + ^-Xl) 

+fe^^ aP - (2.20) 

where we use X\ X a , ip a/ 3, ip a/3 to denote the fluctuations. 
The lowest excited states around the vacuum are 



i 



Tr( ^2^ X/ ~'V^ n/)|0) ' Tn/;W3|0) ' Tr ^ |0) - (2 ' 21) 

The bosonic excitations have energy 

E = \%, (2.22) 

which depends only on the background geometry. And the fermionic excita- 
tions have energy 

E = fe+, (2.23) 

which depend on the background fluxes. If <fi = ax + + b, a, b are constant, the 
first excited bosonic and fermionic states will have the same energy. These 
energy are generically time- dependent except a few special cases. One of such 
exception is the usual plane-wave background, which has been discussed in 
|llj . Another special example is 

= -cr, / = f e CT , (2.24) 
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with X — fo— const., so the first excited states have the fixed energy. 

The higher excited states can be discussed similarly. As the first excited 
states, their energy will be generically time-dependent. 

Though we have computed the energy of the excited states, they are only 
perturbational results in the matrix model. They are only trustable if the 
coupling in the matrix model be very small. If we define 



x^J—x 1 , n / = v /^Tn / , (2.25) 




the bosonic part of the matrix model Hamiltonian about this vacuum is 

Htos = Tr{^fl? + \x] 

+ (R-e-i ) 4 [X / , X J , X K , £ 5 ] [X 1 , X J , X K , C 5 ] (2.26) 
-^{R_f^) 2 {e ijkl X l [X j , X\ X\ £ 5 ] + e abcd X a [X b , X c , X d , C 5 ])} 
The 't Hooft coupling in fl2T26l) is 

g = JRle-f = — (2.27) 

In the early time, r — > —00, the 't Hooft coupling will be very weak. So 
we can trust the perturbation theory of matrix model. On the contrary, in 
the late time, when r — > 00, the perturbation matrix calculation will be not 
convincing, but the perturbative string theory is well-defined then. 
For the case (12.261) . defining 



' f e -cr / R_ 



°-V&'- ! \/^ n " < 2 ' 28 > 

we obtain the normal-ordered free Hamiltonian 

H® = Tr[/e- CT (4a, + + ^^)\. (2.29) 

The zero-point energy is canceled between bosonic and fermionic sectors in 
the background ( 11.11) . The free Hamiltonian is the same as the one with 
only zero modes in perturbative string vacuum. This suggests that in the 
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string theory limit, the static vacuum becomes the vacuum of string theory 
on plane-wave background with linear null dilaton. 

In the tiny graviton matrix model on plane- wave background, if one takes 
the string theory limit, the 't Hooft coupling becomes very large, which 
suggests that the closed string will appear as non-perturbative objects of 
the model. In our case, the string theory limit is more subtle. Since the 
background is time-dependent, the matrix model and its 't Hooft coupling is 
also time-dependent. Strictly speaking, the matrix model is only well-defined 
in the very early time when the matrix coupling is small and the string 
coupling is large, and on the other hand the string theory is well-defined in 
the very late time when the string coupling is small and the matrix coupling 
is large. If we take the string theory limit at fixed time where the string 
coupling is not very large, the 't Hooft coupling becomes very large just like 
the plane-wave background cases, and the vacuum in matrix model changes 
to closed string vacuum.. 

3 The fuzzy three sphere solutions 

In this section, we will discuss fuzzy three sphere solutions. In our time- 
dependent matrix model, there is no static fuzzy sphere solution and the 
only possible fuzzy three sphere solutions are time-dependent. Since our 
background keeps SO (A) x SO (A) symmetry, without losing generality, we 
take the following ansatz 



x a = ^ = = o, 

X* = S(r)X l , 



i = 1,2,3,4, a = 5,6,7,8. 




Here Xq 7 s are constant J x J matrices satisfying 



[Xq,X j 0: Xq,£ 5 ] 




i 





i=l 



such that 



[X\X\X\C,] 



S 2 (r)e ijkl X l , 



i 




(3.3) 



i=l 
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With this ansatz, the equation of motion for X % is 

4R f 

X i + XX i ^e ijkl [X j , x k , x\ c 5 ] 

f>2 -2<p 

^— [[X\X\X k ,C 5 },X\X k ,C 5 ] = 0, (3.4) 

which gives the equation for S(t), 

S(t) + XS- AR^fS 3 (r) + 3i£e~ 2 ^ 5 (r) = 0. (3.5) 

The linear term of the above equation stems from the mass-squared term in 
the matrix model, which is determined by the background geometry. The 
other two terms, cubic one and quintic one, are from the interaction terms in 
the matrix model. There are two parameters 4>(t) and /(r) in the equation. 
In this article, we discuss two cases whose corresponding perturbative string 
have been studied in [T9] : 

• Case 1:0= — cr, /(r) = foe CT , c > 0, f = const; 

• Case 2: = — cr, /(r) = f , c> 0, f = const. 
In both cases, = — cr, so (13.51) is 

S(t) + f 2 e~ 2cT S - AR_fS 3 (r) + 3R 2 _e 2cT S 5 (t) = 0. (3.6) 

It is very difficult to obtain analytic solutions of this equation. Using the 
numerical method to discuss the solution and combining the analysis of the 
early time and late time asymptotic behavior of the solution, we will obtain 
a good qualitative understanding of the solution. 

For the first case, A = /q is a constant. The equation for S(t) is 

S(t) + / 2 5(r) - 4i?_/ e cr S 3 (r) + 3Rl e 2cr 5 5 (r) = 0. (3.7) 

The solutions depend on three parameters, R-,fo,c. For a definite discus- 
sion, at first, we let 

i?_ = 2, / = 1, c=l. (3.8) 

The figure [1] is for the numerical solution between r = —5 and r = 5 if 
we set initial conditions 

S(0) = 1, 5(0) = 0. (3.9) 
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-1.0-| , 1 , 1 , 1 , , , , , 1 

-6 -4 -2 2 4 6 



Figure 1: The solution of Case 1 when r G [—5, 5] and i?_ = 2, f — 1, c — 1. 
At late time, this solution vibrate faster with a smaller amplitude. In fact, 
its amplitude will shrink to zero very quickly. 

From this figure, we can read out the behavior of S(t) in [—5,5]. When r 
becomes larger, in other words, when g s becomes weaker, S(t) vibrates more 
rapidly, and the amplitude becomes smaller. The figure can not tell us the 
whole story, because when r 3> 1, S(r) vibrates so rapidly that the computer 
can not give right numerical solutions. Fortunately, when |r| is large, we can 
analyze the solutions from their asymptotic behaviors. 

At early time, if we take r — > — oo limit, the equation is dominated by 
the linear term and is approximated to be 

S(r) + / 2 S(r) = 0, (3.10) 

which has the solution 

S(r) = Asm(f r + B). (3.11) 

A, B are constants which depend on initial conditions. So at early time, fuzzy 
three sphere is in a stable vibrating state. 

At late time, when r — > oo, the equation (13. 7p is dominated by the 
interaction term and can be approximated to be 

S(t) + 3i?^e 2cr ,S 5 (T) = 0. (3.12) 

So at very late time, S(t) will vibrate very rapidly around zero. The am- 
plitude will also decay to zero very fast, at the rate of e~^ T . This behavior 
is expected. From matrix model point of view, at late time, the coupling 
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is stronger and the nonperturbative effect becomes essential and the funda- 
mental closed string degree of freedom becomes important so the fuzzy effect 
is quite weak. 

Therefore, we have the following picture. A fuzzy three sphere, vibrates 
in a stable rhythm in the early time. But after a long time evolution, at late 
time, the non-Abelian degrees of freedom of the matrix model become less 
important and the fundamental string degrees of freedom take it over. 

Changing parameters in the equation will not change the whole picture 
but can postpone or accelerate the solutions to reach their asymptotic be- 
haviors. If we change the setting, let 

R- = 10, / = 1, c=l, (3.13) 

and still set the initial conditions 

S{0) = 1, 5(0) = 0, (3.14) 

we obtain the figure El 



2 - 



-1 - 




-6 -4 -2 2 4 6 



Figure 2: The solution of Case 1 when r G [—5, 5] and i2_ = 10, f = 1, c = 1. 
At late time, it vibrates more rapidly than the one in i?_ = 2 case which is 
plotted in figure [U 



Comparing the figures ([!}) and ([2]), we find that at late time, S(r) vibrates 
more rapidly when i?_ becomes larger. This is because the coefficient of 
the quintic term is larger, which suggests a stronger interaction and leads to 
faster vibration. In other words, a larger i?_ makes fuzzy effect weaker. If 
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we let fo be larger, from (13.71) . we can deduce that the fuzzy three sphere 
will shrink more slowly. 

For the second case, A = foe~ 2cT , the equation ( 13. 5ft is 

S(r) + fle~ 2cT S - AR.f S 3 {T) + 3R 2 _e 2cT S 5 (r) = 0. (3.15) 

At first, we set i?_ = 2, / = 1, c = 1 just like (I3.8p . If we input 

5(0) = 1,5(0) = 0, (3.16) 

we have the figure [3] for r £ [—5,5]. 




Figure 3: The solution of Case 2 when r £ [—5, 5] and R- — 2, f — 1, c = 1. 
At late time, its asymptotic behavior is similar to the asymptotic behavior 
of solutions in Case 1. At early time, the solution has a similar asymptotic 
behavior to its late time asymptotic behavior. 



From this figure, we find that the late time asymptotic behavior of this 
solution is very similar to the one in Case 1. But the early time asymptotic 
behavior is very different. Just like Case 1, we can not obtain conclusions 
from figure when |r| 3> 1. So we must analyze the asymptotic behaviors of 
the solutions. 

At late time, when r — > oo, the equation (13.151) is approximated to be 

S{t) + 3R 2 _ e 2cT S\r) = 0. (3.17) 

It is the same equation as (I3.12p . So the asymptotic behaviors of solutions in 
these two cases at late time are the same. We know at late time, A, g s — > 0, 
f = fo still is constant. So the decay of fuzzy three sphere is a natural phe- 
nomenon because the background now is nearly flat and the string coupling 
is weak and the matrix coupling is strong. 
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At very early time, the equation (13. 15j) is approximated to be 



S(r) + fle~ 2T S(T) = 0. (3.18) 

This equation tell us that when r — > — oo, the asymptotic behavior of the 
solution of ( 13.151) will be very similar to its late time one. So the behav- 
ior shown in figure ([3]) is inevitable. To understand these behavior, let us 
come back to equation ( 13. 4ft . The most important term dominating the 
early time behavior, is XX 1 , which depends on the background geometry. 
Since lim r _ > „ 00 A = oo in this case, the large spatial extension reduce fuzzy 
effect. This phenomenon happens also in the matrix models on other time- 
dependent backgrounds [51 124]. 

Nevertheless, we can read out the whole behavior of solutions from the 
figure [31 (13.171) and (13.181) . If there is a very small fuzzy three sphere at early 
time, it will grow up as time goes by. And after it reaches its maximum, it 
becomes smaller and smaller. At late time, it vibrates fiercely and decays 
out. 

We can also analyze the effect of changing parameters. If we use the 
setting 

A_ = 10, / = 1, c=l, (3.19) 

we get the figure H] for the solution between r = —5 and r = 5 with the initial 
conditions (I3.16|) . 



As before, a larger R_ makes the solution vibrate more rapidly when r — > oo, 
and makes the solution vibrate more slowly when r — ► — oo (Since we only 
plot the figure when r G [—5,5], the early time asymptotic behavior of the 
solution is not evident in the figure 4). If we let fo take a larger value, 
however, we will find contrary effect. 

Before ending our study of fuzzy three sphere, it would be interesting to 
discuss how the tiny graviton scale changes with time. In the appendix, we 
investigate the giant gravitons in the background ( 11.11) . We learn that the 
size of giant graviton is changing with time: 

^ + Ar - 4/( -f)r 3 + 3(^)Ve-^ = 0, (3.20) 

where \x = 27r 2 T 3 with T 3 being the tension of D3-brane and P_ is conserved 
lightcone momentum. In the case <fi = constant, A = /g, there exist static 
giant gravitons with radius r\ = fo — and Tq = 3/q — . The former one is the 
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-10-1 , 1 1 1 , 1 , 1 1 1 , 1 

-6 -4 -2 2 4 6 



Figure 4: The solution of Case 2 when r G [—5, 5] and i2_ = 10, fo = 1, c = 1. 
At late time, this solution shrinks to zero much faster than the one when 
we set R_ = 2,f = l,c = 1. At early time, though its full asymptotic 
behavior can not be seen from the figure because we do not plot the figure 
of the solution when r < —5, we can confirm that its asymptotic behavior is 
similar to i?_ = 2 case through analysis in this section. 



well-studied stable giant graviton keeping half supersymmetries. The latter 
one is not stable. From this, one can read out the size of tiny graviton which 
carry one unit of lightcone momentum, 




' tiny • \ <J-i '- 1 -/ 

^ AdS 

For the case being not a constant, after identifying P_ = fi/R-, we can 
find that the equation (13.201) takes the same form as the equation (13.51) . This 
shows that the scale of tiny graviton changes with r in the same way as the 
size of fuzzy three sphere. This is not a surprise since the tiny graviton is the 
fundamental parton of our matrix model. This fact could be understood as 
following: the fuzzy three sphere could be taken as the blow-up of the tiny 
gravitons, therefore its size is changing with the evolution of the scale of tiny 
graviton. 

Actually the fact that the dynamics of fuzzy three sphere is in perfect 
agreement with the dynamics of giant graviton in the background fll.il) gives 
a consistent check that our matrix model on time-dependent background is 
correct. 
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4 Conclusions and discussions 



In this paper, we constructed a matrix model for type IIB superstring 
theory on a plane- wave background with null dilaton (11. ip . Our matrix model 
is time-dependent, have much different character from the tiny matrix model 
on plane- wave background. One of the differences is that the vacuum solution 
is the unique static solution. The spectrum about this vacuum are time- 
dependent in most generic cases. We compared our second order Hamiltonian 
of trivial vacuum to Hamiltonian of string er-model on the linear null dilaton 
background, and found that the spectrum of perturbation of background in 
string theory could be produced from our second order Hamiltonian. 

We also investigated the fuzzy three-sphere solutions in our matrix model 
and found that they must be time-dependent. The solutions depend on the 
background geometry and flux. In this paper, we discuss two cases with 
linear null dilaton. One is = — cr, /(r) = foe CT , f = constant. The other 
one is (f> = — cr, /(r) = fo = constant. For these two cases, we could not 
obtain analytic solutions yet, but we investigated the solutions numerically 
and also analyzed their early and late time asymptotic behaviors. The results 
are interesting. In the first fuzzy sphere vibrates stably at early time, 

but at late time it vibrates much more rapidly, then decays out. In the second 
case, at early time a very small fast- vibrating fuzzy sphere appears and grows 
up gradually to reach its steady maximum and then decays through rapid 
vibration at late time. The picture can be understood as following. The early 
time behavior of the fuzzy three-sphere is determined by the background 
geometry completely, so in the second case, the fuzzy degrees of freedom is 
diluted. On the other hand, the late time asymptotic behavior is determined 
by the interaction terms in the matrix model. The very strong interaction at 
the late time requires the nonperturbative degrees of freedom to replace the 
matrix degrees of freedom. This explain why the fuzzy sphere always decays 
out at the late time. 

The matrix model constructed here is well-defined at the very early time 
where the coupling of the interaction is weak. It is dual to the late time 
perturbative string description. In other words, the closed string perturbative 
string will appear as the nonperturbative objects in the matrix model, and 
vice versa. The dual pictures give a complete description of the string theory 
in the linear null dilaton background. 

In the matrix model on plane- wave background [TT], the spectrum around 
fuzzy sphere vacuum contain photon states, which are the gauge field on the 
giant graviton in the string theory limit. In our case, fuzzy sphere solutions 
are too complex and time-dependent. We did not go to similar discussions 
in this paper. We hope we can discuss these issues in the future. 
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One essential ingredient in the construction of tiny graviton matrix model 
is the existence of Nambu three-bracket in the D3-brane action, which allow 
us to do matrix regularization to obtain a matrix model. On the other hand, 
the Nambu three-bracket is a realization of infinite dimensional 3-algebra, 
which plays a central role in the recent study of multiple membrane theory 
by Bagger-Lambert [2S] and Gustavsson[2"T]. In [2E1 I2H], it was pointed 
out that the space of functions on three dimensional manifold, with Nambu 
three bracket could be identified as an infinite dimensional 3-algebra. The 
Nambu three bracket is the realization of three-bracket of this infinite di- 
mensional 3-algebra. Especially, it was proposed in [30J that the deformed 
BLG-theory on R x T 2 gives tiny graviton matrix theory of Type IIB plane- 
wavej^] Very recently, in [32], the authors used the similar way to [TT] as an 
matrix realization of "relaxed three-algebra ", which is a '3-algebra' with re- 
laxed closure and fundamental identity conditions. The relaxed three-algebra 
is a linear space of matrices with a relaxed three bracket, which is defined 
as [T a ,T b ,T c ,T~] with T~ being a special fixed constant matrix. Using re- 
laxed three-algebra, they tried to construct a Hermitian model to describe n 
(which is a general number) membranes systems. These discoveries suggest a 
deep but not clear relation between multi-membrane theory and tiny graviton 
matrix model. It would be important to understand this issue better. And 
it is also interesting to discuss the multiple membrane in a time-dependent 
background [33]. 
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Appendix: Giant gravitons 

In this appendix, we would like to study the giant gravitons in the back- 
ground ( 11. ip . Let us introduce the polar coordinates (r, 9, </>, cp) in the plane 
defined by the Cartesian coordinates Xi,i = 1,2,3,4. In polar coordinates, 
the 5-form flux 

(F 5 ) +12U = (F 5 ) +r ^ = 4/(x + )r 3 sin 2 9 sin (4.1) 

3 The relation between mass deformed M2 theory on T 2 and IIB DLCQ string on plane- 
wave was first pointed out and identified in |31j . 
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with the relevant nonvanishing components of 4-form potential 



C+^ = /(x + )r 4 sin 2 £sin0. (4.2) 

Here we just focus on the case that the giant graviton being blown up in 
Xi,i = 1,2,3,4. Obviously the case that the giant graviton being blown up 
in Xi, i = 5, 6, 7, 8 could be studied similarly. 

The giant gravitons correspond to spherical D3-brane embedded in the 
background consistently. The D3-brane is wrapped around S 3 characterized 
by 9, 0, if. This allows us to identify the worldvolume coordinates of D3- 
branes to be r, 9, <p and </?. To simplify the discussion, we can choose the 
gauge X + = t such that the only nontrivial embedding is r depending on r. 
After this truncation, the action of the D3-brane reduce to that of a point 
particle: 

S = V J dr[-e-^r 3 (r) \/ -G AB d T X*d T X* + /(r)r 4 (r)], (4.3) 

where /i = 2n 2 T 3 E As we have seen, the translation alone; x is still an 
isometry so its corresponding canonical momentum is a conserved quantity. 
And the canonical Hamiltonian is just H = —P+, which is of the form 

2P_ 2P_ 2 p K J 

where P r is the canonical momentum of r. The Hamiltonian describes a par- 
ticle with mass P_ moving in a time-dependent potential. The r is actually 
the size of the giant graviton, which is determined by the equation of motion 

tl- + Ar - 4/^r 3 + 3(-^)V 2 V = 0. (4.5) 

CLT Z P- P_ 

This equation is the same as (13. 5p after identification R_ = fi/P-. Therefore 
the dynamics of giant graviton is in exact agreement with fuzzy three sphere. 
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